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We study the following algebraic system of equations 

Pa(xi, . . . ,Xn) := Paix) := x"" -^ac,/3X^ = 0, a e J, (1) 

where I = {a E Z"^ : \a\ < m}, J := Jm+i = {« G Z!^ : |a| = m + 1}. 
We consider some i^I x ^/ - matrices: ^41, . . . , A^. The structure of these matrices is as 
follows: k rows {k = i^Jm) are the coefficients of k polynomials from {Pa ct E J} and 
each of the remaining rows consists of O's and one 1. 

We prove that the number of distinct solutions to (1) equals to maximal possible if 
and only if the above matrices are commuting and each of them is semisimple (i.e. they 
have complete system of eigenvectors). 

This gives a characterization of correct systems of points for multivariate Lagrange 
interpolation. This also reduces the solving of (1) to solving of n univariate algebraic equa- 
tions: to finding the eigenvalues of the above matrices. In the case when the eigenvalues 
of some matrix are distinct, the solving of (1) is reducing to solving of only one univariate 
equation: to finding of these distinct eigenvalues. 

The above results remain valid also for more general systems (1): with / being any 
lower set. 


